March 2005 



ON THE COHERENT STATE PATH INTEGRAL 
FOR LINEAR SYSTEMS 

C. G. Torre 

Department of Physics, Utah State University, Logan, Utah, 84322-4415 USA 

in 

' We present a computation of the coherent state path integral for a genèric linear system using "functional 

çSJ , methods" (as opposed to discrete time approaches). The Gaussian phase space path integral is formally given 

5^ ' by a determinant built from a first-order differential operator with coherent state boundary conditions. 

We show how this determinant can be expressed in terms of the symplectic transformation generated by 
the (in general, time-dependent) quadratic Hamiltonian for the system. We briefly discuss the conditions 
under which the coherent state path integral for a linear system actually exists. A necessary — but not 
| sufhcient — condition for existence of the path integral is that the symplectic transformation generated by 

the Hamiltonian is (unitarily) implementable on the Fock space for the system. 
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> 
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1. Introduction 

The coherent state path integral, a variant of the phase space path integral, has long 

m 
o 



been recognized as a useful tool in quantum mechanics and in quantum field theory (see, 
e.g., [l]-[5]). In the quantum mechanical setting this form of the path integral has been 
studied fairly extensively (see, e.g., [4]), although it would seem that only relatively re- 



cently does there appear a definitive treatment of the coherent state path integral for a 
genèric linear system with one degree of freedom (in the context of a study of the semi- 
classical approximation) [5]. In quantum mechanics the method normally used to define 
and analyze the coherent state path integral is based upon taking a limit of an approx- 
imation based upon discretized paths. In the field theoretic setting the coherent state 
path integral is mainly used to set up the perturbative evaluation of the S-matrix; here 
"functional methods" (as opposed to discretized path methods) are normally employed. 
In particular, at the level of the free or semi-classical theory the path integral represen- 
tation of the vacuum to vacuum transition amplitude, after integrating out the canonical 
momenta, is expressed in terms of a Fredholm determinant of a linear differential operator 
(wave operator, Dirac operator, etc. ), which features in the approximate quadratic action 
functional. 

In this paper we compute the coherent state path integral for a genèric linear bosonic 
dynamical system with general coherent state boundary conditions. The number of degrees 
of freedom can be infinite, so this computation includes field theory. In our computation 
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we do not first integrate out the canonical momenta, nor do we restrict attention to rela- 
tivistic fields. Therefore, the path integral is formally given in terms of the determinant 
of a first-order differential operator with coherent state boundary conditions. We show 
how this determinant can be expressed in terms of the symplectic transformation gener- 
ated by the (in general, time dependent) quadratic Hamiltonian for the system. While 
our computations are somewhat formal, the resulting expression for the coherent state 
transition amplitude agrees with the (rigorous) result obtained using methods of canonical 
quantization [6], as it should. 

The computation provided here demonstrates a viable method for evaluating a class 
of coherent state path integrals which we hope will be a useful addition to the collection 
of techniques used to compute path integrals. This computation should be relevant for a 
number of applications, including: any linear quantum mechanical system, linearized/semi- 
classical approximations, quantum fields in curved spacetime, and parametrized free field 
theory and various other quantum gravity models [7]. 

Our results are also intended to help explore within the path integral formalism subtle 
quantum field theoretic phenomena that have been uncovered using other methods of 
quantization (e.g., canonical quantization). In particular, it is well-known that for field 
theories there exist inequivalent representations of the canonical commutation relations 
(see, e.g., [8] and references therein). Closely related to this is the fact that many linear 
canonical transformations — which may include those defining the time evolution of a linear 
system — cannot be unitarily implemented in the Fock space quantization of a field theory 
[9,10,6]. This situation is known to occur in a variety of physical settings, e.g., for quantum 
fields in curved spacetimes [11], for the polarized Gowdy model in General Relativity [12], 
and for parametrized free field theories in dimensions greater than two [13]. It is natural 
to ask how these important phenomena, which have been understood heretofore using 
operator techniques, manifest themselves in the path integral formalism. Because we can 
express the path integral in terms of the symplectic transformations generated by the 
classical Hamiltonian, the connection with results from canonical quantization on unitary 
implementability/equivalence becomes accessible. For example, in this paper we shall 
see explicitly that unitary implementability of dynamical evolution is necessary but not 
sufficient for the coherent state path integral to exist. 

2. Preliminaries 

We will be considering the path integral for a linear dynamical system, by which we 
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mean the following. Fix a real Hilbert space, i.e., a real vector space V, complete with 
respect to a scalar product (•,•). Elements of V will be denoted z, w, etc. The vector 
space V is to be the phase space for the system, so we further assume that V is equipped 
with a densely defined symplectic form O and a Hamiltonian 7i, which is a densely defined 
quadratic form on V. As explained e.g., in [14], we require the inner product and symplectic 
form satisfy 

v ' 1 4 w ^ (w, w) K 1 

This implies that the symplectic form is bounded on V and so can be defined with V as 
its domain. The scalar product and symplectic form then combine to define a (bounded, 
skew-adjoint) complex structure J: V — > V via 

Q(z,w) = 2(z,Jw), (2.2) 

which can be used to define "positive and negative frequency" solutions to the linear field 
equations defined by Tí. More precisely, we introduce the complexification of V and 
extend (•,•), O and J to V c via linearity. We define complex conjugate Hilbert spaces 
V± C corresponding to the ±i eigenspaces of J equipped with the sesquilinear inner 
products 

z,w G V± ->2(z *,«;). (2.3) 

Here we use an asterisk for complex conjugation on . For later convenience, we suppose 
we have fixed a basis: {e a } for V+ and {e^} for V-. If z E V, we have* 

£=^e a + ^*e*, (2.4) 

a relationship we will often denote as either 

( 2 - 5 ) 

or simply 

If z G we write 

z = z a e a + z a e*, (2.7) 



* The summation convention is in effect and the range of the summation can be infinite and 
continuous (i.e., integration) . 



or 



or simply 



\a ) , (2-8) 



(2-9) 



we have 



z 



z a 

so that z = ( _ a ) is in fact an element of V if and only if z a = z a * . With this notation 



zeV+ < — ► z = z a e a < — ► z=\^ Z J, (2.10) 

zeV- <— > * = <— > (2.11) 

z*eV- <— z = ^*e* <— f=^2^. (2.12) 
The basis {e a , e *} is chosen such that 

(w,z) = ^(w a z a + w a z a ), (2.13) 

Q(w,z) = -(z a w a -w a z a ), (2.14) 
i 

and 

= A flò ^ 6 + ifl^z 6 + is^z 6 . (2.16) 

where A is self-adjoint and B, B = B* are symmetric. 

As explained, e.g., in [14] the data V, O, /j, define a "1-particle" Hilbert space and 
a corresponding Fock space representation of the Heisenberg group. In this representa- 
tion, the operators corresponding to the classical variables z b and z b are annihilation and 
creation operators a b and a'^, respectively. 

The path integral we will be studying is over a space of paths z = z(t) in V . The 
paths will be required to obey coherent state boundary conditions of the form 

z a (0) = w a , z a (T) = v a , (2.17) 

for some given w a and v a . Given the Hamiltonian (2.16), symplectic structure (2.14), and 
boundary conditions (2.17) the action functional on this space of paths is given by 

T dt { ^(z a z a - z a z a ) - (A ab z a z b + L -B ah z a z h + L -B ah z a z h ) \ 

(2.18) 



Q(z) = J dt \^Cz a z a - z a z a ) - (A ab z a z b + \B ab z a z h + l -B ab z a z b )} 



+ ±[z a (T)z a (T)+z a (0)z a (0)]. 



This is just the phase space action functional on V expressed in a complex coordinate chart 
and extended to V C . The boundary terms in (2.18) are there so that Q is differentiable 
with the boundary conditions (2.17) [2]. 

In light of (2.17), we fix a path z(í) with these boundary conditions and set 

z(t) = z(t) + C(t). (2.19) 
The variables (f(í) satisfy the boundary conditions 

C a (0)=0, C(T) = : (2.20) 
and are elements of a vector space X with dual X' . The pairing between X and X' is 

ex, C'(í)= (|) el' - (C\Ò = £ dt^mC a (t) + C(t)C(t)). 

(2.21) 

The action restricted to X is a quadratic form characterized by a symmetric linear operator 

Q(0 = {D(,0, (2.22) 

where 

of= (i|-Y'-y'). (2.2J) 

V "«Ca - A aòC - B abQ / 

We also define the quadratic forms Qo and V on X and symmetric operators Dq: X — > X' 
and N:X^X' by 



Qo(0 = ^ ^{^(C fl C a -CaC a )-44C a 
= <A)C,C>, 



(2.24) 



and 



where 



Q(Ò = Qo(C) + nO, (2-26) 



3. The Path Integral 

Our goal is to compute a path integral which is, roughly, of the form 

e^, (3.1) 



where the integral is over a space of phase space paths with the boundary conditions (2.17). 
From the point of view of the Fock space representation of the quantum theory, I(v, w) 
should be a coherent state matrix element of the time evolution operator £7(0, T) defined 
by a Hamiltonian corresponding to the classical observable H: 

l(v,w) = {v\U(0,T)\w). (3.2) 

Here the states \v) and \w) are Fock space eigenvectors of the annihilation operator a b : 

a b \v} = v b \v), a b \w) = w b \w). (3.3) 

An approach to obtaining a rigorous definition of (3.1) can be found in [15,16]. We 
shall not attempt to define functional integration here; instead we shall proceed formally 
by postulating a few bàsic properties, which feature in the rigorous definitions of [15,16], 
and which any suitable definition of integration should exhibit. We shall assume that 
integration is a linear operation on a class of functions on the vector space X. The 
integration operation is normalized relative to the quadratic form Qq and is denoted with 
the symbol Jx[dÇ]Q - The normalization condition is 



= 1- (3-4) 

X 



JQo(C) = 

We formally define the path integral (3.1) by 

l(v,w)= í [dC] Qo é Q ^. (3.5) 
J X 

As we shall see, in the Fock space representation our choice of path integral normalization 
corresponds to defining the time evolution operator using the normal-ordered Hamiltonian 
operator associated to the classical expression (2.16), with no additive c-number renormal- 
izations. 



To evaluate the path integral we choose 



z(t) 



^ = \W)) (3 - 6) 

to be a critical point of the action functional, i.e., z(t) is chosen to be the (unique) path 
in satisfying 

Dz(t) = 0, z(0) = w, z(T) = v. (3.7) 
Using linearity of the integration operation, the path integral now takes the form 

l(v,w) = exp j^ a z a (T) +z fl (0K]} J x [dC] Qo e iQ ^ . (3.8) 
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The remaining path integral is an oscillating Gaussian — or Fresnel - type of integral. It 
is (at least formally) given in terms of the Fredholm determinant [16]:* 

J(0, 0) = í [d(] Qo e^O = det-kD^D). (3.9) 
J X 



Hence 



X(v, w) = exp |^Kz a (T) + z a (0K] j det^L»- 1 !)). (3.10) 



4. Evaluating the path integral in terms of classical dynamics 

We now show that the result (3.10) can be expressed in terms of the symplectic trans- 
formation(s) generated by H(t) and the boundary conditions (2.17). We begin by giving 
some relevant features of the symplectic transformation. 

We assume the Hamiltonian H(t), which is allowed to be time-dependent, generates a 
1-parameter family of symplectic transformations S(t) on . S is determined by 

± + H(t))S(t) = t H(t) = i(*$ t) _^)) ; 5(0) = id. (4.1) 

The operators S(t) can be expressed as a Dyson-type of expansion, i.e., the time-ordered 
exponential of the linear transformations H(t) on V C defined by the Hamiltonian H(t). 
We shall not need the explicit formula here. Using a block matrix notation paralleling 
(2.8), we can express the symplectic transformations as 

i.e., 

— I 

where the Bogoliubov coefficients a(t), (3{t) satisfy at each t 

a:V+^V+, /3:V-^V+, (4.4) 

a:V-->V_, P:V + ^V-, (4.5) 
a = a\ P = P*, (4.6) 



* In principle there could also be a phase factor coming from the index of Q relative to Qq 
[16], but we shall see that this factor is unity. 
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K t _ = i( i v+j a (3 T - (3 a T = o. (4.7) 
From these equations it follows that a(t) has a bounded inverse and that 

Finally, we have that 

a(0) = id v+ , a(0) = zrf v _, /3(0) = 0, (3(0) = 0. (4.9) 

All the dependence of I(v,w) on the initial and final states is in the exponential of 
the action sitting in front of the determinant in (3.10). To make this initial/final state 
dependence explicit, we use the fact that the critical points 

«•«» 

of the action are determined by the symplectic transformations generated by the Hamilto- 
nian (2.16): 

z(£) = «S(í)z(0). (4.11) 
The boundary conditions (2.17) imply 

z(T) = a~^(T)w + a(T)v (4.12) 

z(0) = a _1 (T)t; -^(T)w, (4.13) 
where we have defined symmetric operators 7: V+ — > V- and a: V- — > V+ via 

7 = a- 1 / 3, a = /3a- 1 . (4.14) 

Putting this all together, we get 
J(t/, «;) = exp {(a- 1 ï) a b (T)v a w b + l -a a b (T)v a v b - ^ 7 J(7> a «; 6 } ^{D^D). (4.15) 

The determinant appearing in (4.15) depends upon the Hamiltonian (2.16) via the 
symplectic transformation S. We will now make this dependence on S explicit. The 
hypothesis on M = Dq 1 D is that M — 1 is trace class; this guarantees the determinant is 
well defined and satisfies the variational identity. 

ílogdet(M) = Tr(M _1 <5M). (4.16) 
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The trace ( "Tr" ) of an operator R: X — > X is defined in terms of an integral kernel R(t, u) 
with vàlues in the set of operators on V C : 

Rz(t)= [ duR(t,u)z(u). (4.17) 
Jo 

We set 

Tr(R) = I dttvR(t,t), (4.18) 
JO 

where "tr" denotes the Hilbert space trace. To use (4.16) we define the differential operator 
Dy.X -> X' by 

D X Ç=( ^-A ba C b -XB ab ? \ 

(so that D = Di). Setting 

J A = det-^( J D c 7 1 J D A ), (4.20) 

we have 

^logX A = -^{{D^D^D^n} = -^(D^N), (4.21) 

where N was defined in (2.27). Our strategy is to obtain a suitable expression ofTr (£>-%) 
in terms of symplectic transformations and then solve the differential equation (4.21) in À 
with initial condition given by the path integral normalization, 

J = 1, (4.22) 

to find X A , from which we have X = I\ = i* 

To obtain a suitable form of T^D^iV) we need to obtain an expression for the Green 
function G\ = D^, which is uniquely determined by the solution of the system 

D x Ç(t)=F(t), C a (0)=0, C(T) = 0, (4.23) 



so that 

rT 



C(t) = G\F(t) = ! du G x (t, u)F(u). (4.24) 
Jo 

The general solution to the differential equation in (4.23) is easily checked to be 

C(í) = Co(í)+/ du6(t-u)S x (t)S^(u)i:F(u), (4.25) 
Jo 



* This approach will define X provided D\ has no zero eigenvalues as À varies from to 1 
[16]. Otherwise there is an additional phase coming from the index of D\. We shall see 
that D\ has no kernel. 



— * _ —* 

where Çq is the general solution to -D^Co 

Co(í) 



= o, 

= <s a (í)Co(o), 



(4.26) 



with 



(4.27) 



being the symplectic transformation generated by 



H x (t) = A ab (t)z a z b + i\ (B ab (t)z a z b + B ab (t)z a z b ) . (4.28) 



S\(t) satisfies 



(4.29) 



and, for each value of À, the obvious generalizations of equations (4.4)-(4.7) and (4.9) hold 
In (4.25) 9 is the step-function 

9(x) = { 

and 

' -i 
z 



1 x > 0, 
x < 0, 



The boundary conditions in (4.23) imply (suppressing indices) 



where 



Co(f) 



Co(0) = -a^(T) 



( &(*)Co(0) , 

v«a(*)co(o);' 



rf-u l S A (T)5^ 1 («)SF(-u) 



Here we use a notation for the components of elements of such that if 



then 



(v) ( = V, (v)ç = V. 

Putting all of this together, the solution 

Mi) 
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(4.30) 
(4.31) 



(4.32) 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



to (4.23) is given byf 



C(t) =J o du{ (0(t - u)a x (t) - /3 A (í) 7A (T)) 
-6(u-t)(3 x (t) [5^(w)SF(«) 

C(t) =J o du{ (e(t - u)P x (t) - a x (t) lx (T)) 
-6{u-t)a x {t) Í5~ 1 (m)EF(w) 



S^(u)EF(u) 



(4.37) 



where 



C 



-ia[F c - iplFç, 



S^EF 



ií3{F c + ia T x F- c . 



The Green function thus takes the form 



(4.38) 



(4.39) 



where 



-i9(t-u){a x (t)a\(u)-p x (t)pl(u)) + i(3 x (t) (l X (T)a\(u) - p{(u)) , 



iO(t - u) ( a x {t)(3 T x {u) - (3 x (t)a{ (u) ) + i(3 x {t)[ lx {T)(3{ («) - a{ (u 



T 



T i 



G^(t,u) = -i6{t- u)(p x {t)a\{u) -a x (t)0{(uj) + ia x (t)(~f X (T)a[(u) - p[(uj), 
G rc (t, u) = -iO(t - u) (p x (t)pl(u) - a x (t)a T x (u)) + ia x (t) ( 7 a(T)/3a («) - «a ( 



lu. 



(4.40) 



We then have 



Tr(G x N) = -j dttr{G^(t,t)B(t) + G^(t,t)B(t)} 

= -i /" T rfítr{/3 A (í)( 7A (T)/3f(í) -a T x (t))B(t) + a x (t)(^ x (T)a[(t) - p{(t))É(t) 

(4.41) 



Equation (4.41) can be simplified considerably. In particular, we claim that 



Tr(G A iV)=tr(a A 1 (T)^a A (T)), 



(4.42) 



f As can be seen from this result, there are no non-trivial solutions to -D A C = with boundary 
conditions (2.20). This means that D x has no kernel and there are no additional phases to 
be computed (see [16]). Evidently, this is a simplifying feature of the coherent state path 
integral. 



11 



where the trace on the right hand side of the equation is on V-. To prové (4.42), we vary 
the equations (4.29) determining S x so that 



(4.43) 



Equation (4.43) can be viewed as an inhomogeneous equation for dS x /dX with "source" 
given by — ~^S X . With initial condition 

'dS x 



dX ) t _ 



= 0, 



í=0 



the solution is 
Using (4.27) and 
in (4.45) we get 

d\ 



dS X {t) 
dX 



= -J*duS x (t)S-\u) d ^-S x (u). 



dX 



B 
-B 



(4.44) 



(4.45) 



(4.46) 



t 

i I du 



a x {t)a T x {u) - P x (t)Pi(u) B(u)p x (u) 



T, 



+ 



à x (t)Pl(u)-P x (tW X (u) B(u)a x (u 



(4.47) 



We then have 

tr ( 5A -V)^f) 



rT 

i } du tr 



í 

JO 



+ 



a T x (u)- lx (T)[3 T (u)]B(u)[3 x (u) 
(3[{u) - lx (T)a\u)]B{u)à x {u) 



= Tr { D x lN ) • 



(4.48) 



With this result in hand we can obtain the final form for the path integral. Using 

tr («^Cn^P) = ^[logdet(a A (T))], (4-49) 

and (4.21) we have 



^{lo g X A + ^logdet(a A (T))} = 0, 
so that, using (4.22), we have 



det-^iD^D) = Xi = àet-í{à.Q L {T)ài{T)). 
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— 1/ 



(4.50) 



(4.51) 



Finally, from (4.15) and a\ = a we have 

l(v, w) = exp | (a-^) a b (T)v a w b + \a a b (T)v a v h - \l b a {T)w a w h 

(4.52) 

We note that the coherent state path integral result (4.52) for the transition amplitude 
has been rigorously obtained within the Fock space formalism for a particular class of 
normal-ordered, time-independent, quadratic Hamiltonians [6]. 



^ àet(a Q l (T)a{T))' 



5. Discussion 

Evidently, the path integral as computed in (4.52) makes sense provided the exponential 
factor and the determinant of the operator K = a.Q l (T)a(T) exist. For the Fredholm 
determinant of an operator K to be defined it is necessary and sufficient that K — íd is 
a trace-class operator. (In this section u id n denotes the identity operator on V+.) An 
important necessary condition for det(K) to exist is that the Bogoliubov coefhcient f3(T) 
in S(T) defines a Hilbert-Schmidt operator: 

tr[0t (T)0(T)] < oo. (5.1) 

This can be seen by first noting that ao is a unitary operator, which follows from the easily 
established fact that f3g = 0. Then we have 

|det(ií)| 2 = det(KKÏ) = det(id + 1 PP^a ), (5.2) 

where all operators are evaluated at time T. The operator ao(3(3^aQ is trace-class if and 
only if Pfi* is, i.e., P(T) must be Hilbert-Schmidt. When P(T) is Hilbert-Schmidt it 
follows that the exponential factor in (4.52) exists because each of a _1 (T), cr(T) and 
7(T) is bounded. This necessary condition is highlighted here because it is known that the 
symplectic transformation S(T) is unitarily implementable in the Fock space representation 
defined by fi and J if and only if P(T) is Hilbert-Schmidt [9,6].* Thus the path integral, 
normalized using Qq and interpreted using the Fredholm determinant, fails to exist if 
the symplectic transformation generated by the classical Hamiltonian fails to be unitarily 
implementable in the Fock space representation. 

It should be emphasized, however, that even if the symplectic transformation corre- 
sponding to time evolution from t = to t = T is unitarily implemented, this does not 

* There are are number of interesting situations where this Hilbert-Schmidt condition fails 
(see, e.g., [13]-[17]). 
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guarantee that the path integral exists. Implementability means only that the absolute 
value of the determinant is defined — the phase of the determinant may not be defined. 

In the case where the Hamiltonian is time-independent, some sufhcient conditions for 
the Fredholm determinant of K to exist (and hence for the exponential factor to exist as 
well) can be obtained from the results in [10,6]. For example, a relatively simple sufficient 
condition is that the quadratic form on V+ given by B corresponds (via the scalar product 
on V+) to a Hilbert-Schmidt operator: 



From the point of view of the Fock space formulation of the quantum system, this implies 
that the symplectic transformations S(t) form a strongly continuous group and are repre- 
sented (projectively) as a continuous unitary group generated by a Hamiltonian operator, 
which is unique up to addition of a múltiple of the identity. The Hilbert-Schmidt condi- 
tion on B is then equivalent to requiring that the vacuum state of the Fock representation 
defined by fi and J is in the domain of the Hamiltonian. 
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